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ABSTRACT 

Following the completion of cosmic reionization, the mean-free-path of ionizing photons was set by 
a population of Ly-limit absorbers. As the mean-free-path steadily grew, the intensity of the ionizing 
background also grew, thus lowering the residual neutral fraction of hydrogen in ionization equilibrium 
throughout the diffuse intergalactic medium (IGM). Lya photons provide a sensitive probe for tracing 
the distribution of this residual hydrogen at the end of reionization. Here we calculate the cosmic 
variance among different lines-of-sight in the distribution of the mean Lya optical depths. We find 
fractional variations in the effective post-reionization optical depth that are of order unity on a scale 
of ~ 100 co-moving Mpc, in agreement with observations towards high-redshift quasars. Significant 
contributions to these variations are provided by the cosmic variance in the density contrast on the scale 
of the mean- free-path for ionizing photons, and by fluctuations in the ionizing background induced by 
delayed or enhanced structure formation. Cosmic variance results in a highly asymmetric distribution 
of transmission through the IGM, with fractional fluctuations in Lya transmission that are larger than 
in Ly/3 transmission. 

Subject headings: cosmology: theory - galaxies: formation 



1. INTRODUCTION 

It is thought that the appearance of the first galaxies 
and quasars reionized the cosmic hydrogen in the inter- 
galactic medium (IGM) that was left neutral following 
the recombination era at z ~ 10'^. The process started 
with the emergence of ionized regions around galaxies and 
then groups of galaxies. These H II regions eventually 
overlapped, leaving the IGM highly ionized in low-density 
regions surrounding pockets of residual neutral hydrogen 
in high density regions. The process of reionization was 
inhomogeneous (Furlanetto et al. 2004), and the overlap 
should have occurred at different times in different regions 
(Wyithe & Loeb 2004b). However the identification of 
a Gunn-Peterson trough in the spectra of the most dis- 
tant quasars at redshifts of z ^ 6.3-6.4 (Fan et al. 2004) 
hints that the reionization of cosmic hydrogen was com- 
pleted near z ^ 6, about a billion years after the big bang 
(White et al. 2003; Wyithe & Loeb 2004a; Mesinger & 
Haiman 2004). 

Prior to the overlap epoch, the mean-free-path of ioniz- 
ing photons was limited by the size of individual H II re- 
gions. Following overlap, regions of IGM were exposed to 
ionizing radiation from additional sources. This increase in 
the ionizing intensity ionized hydrogen in denser regions, 
and resulted in an increased mean-free-path which in turn 
led to an even higher ionizing intensity. As the ionizing 
intensity in the post-overlap IGM rose, the optical depth 
to Lya absorption decreased (e.g. Fan et al. 2002). Obser- 
vations of the absorption spectra of high redshift quasars 
(Fan et al. 2004) and potentially gamma-ray burst after- 
glows (Barkana & Loeb 2004; Chen et al. 2005) provide 
clues to the evolution of the ionizing radiation field and of 
the evolution of the ionizing photon mean-free-path. 

The post-overlap evolution in the IGM is expected to 



be inhomogeneous. Indeed values of optical depth to 
Lya absorption measured along different lines of sight at 
5 < z < 6 show fluctuations that are of order unity (Becker 
et al. 2001; White et al. 2003; Djorgovski 2004). The 
Lya optical depth therefore varies from place to place at 
any slice of cosmic time. In addition to small scale vari- 
ations in Lya optical depth generated by structures like 
filaments and mini-halos (the so-called Lya forest), there 
would be large scale inhomogeneities in the evolution of 
the ionization state of the IGM. Fluctuations in the trans- 
mission averaged on a large scale are sensitive to both 
small scale structure along the line of sight (Lidz et al. 
2005), and to cosmic variance on larger scales. Here we 
concentrate on the contribution from cosmic variance in 
the smoothed density field. However the presence of high- 
column density absorbers and small scale structure in the 
density field will add additional fluctuations to those de- 
scribed in this paper. Variations in the average transmis- 
sion would be present on the scale of the mean-free-path 
for ionizing photons (~ 10-100 co-moving Mpc). Indeed 
cosmic variance on scales of 10s of Mpc can be significant 
and has not described by contemporary numerical simula- 
tions of the reionization process (Barkana & Loeb 2004). 
To model cosmic variance properly with numerical simula- 
tions requires resolving low-mass galaxies as well as solv- 
ing the radiative transfer equations inside a very large box 
(hundreds of co-moving Mpc) with a high dynamic range 
of scales. Progress is now being made toward this goal 
(Kohler, Gnedin, Miralda-Escude & Shaver 2005; Kohler, 
Gnedin & Hamilton 2005; Iliev et al. 2005). 

In this paper we examine the variations in the opti- 
cal depth averaged over a redshift bin of sufficient size to 
smooth out the small-scale inhomogeneities of the Lya for- 
est. Our model for the variation in optical depth includes 
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several contributions: (i) variations in the density-contrast 
on the scale of the mean- free-path for ionizing photons, (ii) 
variations in the ionizing background on the same scale 
due to enhanced or delayed structure- formation, (iii) vari- 
ations in the ionizing background due to finite light travel 
time across the clustering length of ionizing sources, (iv) 
variations in the ionizing background due to Poisson fluc- 
tuations in the number of ionizing sources, and (v) evo- 
lution of the density of Ly-limit absorbers. Throughout 
the paper we adopt the set of cosmological parameters de- 
termined by the Wilkinson Microwave Anisotropy Probe 
(WMAP, Spergel et al. 2003), namely mass density pa- 
rameters of ilm = 0.27 in matter, fif, = 0.044 in baryons, 
l^A = 0.73 in a cosmological constant, and a Hubble con- 
stant of Ho = 71 kms~^ Mpc~^ 

2. FLUCTUATIONS IN IONIZING INTENSITY AND LYa 
OPTICAL DEPTH 

The background intensity of ionizing radiation is pro- 
portional to the co-moving luminosity density due to ion- 
izing sources (risic), the density contrast (A), the bias of 
sources in an overdense region (Bgrc), the co- moving mean- 
free-path for ionizing photons (i?mfp), and a factor (Fq) 
describing the fractional contribution of quasars to the ion- 
izing back-ground relative to stars 

J,(Xn,rc(l + ^)'Ai?src^(l + FQ). (1) 

1 + z 

Note that we assume a \mivcrsal spcctriim for the contri- 
bution to arising from a fixed redshift because we are 
averaging over a large number of sources. However the 
contribution to the ionizing background at a fixed point 
in the IGM from redshifted photons will result in a value 
for Ji, that depends on the shape of this universal spec- 
trum, as well as the mean-free-path. On the other hand 
this dependence on contribution from redshifted photons 
is the same in all regions of size i?mfp- Therefore since in 
this paper our interest lies in relative fluctuations of the 
ionizing back-ground among random regions of fixed size 
centered at the same redshift, we assume the exact shape 
of the spectrum to be unimportant for our calculation. 

The mean-free-path (i?mfp) is in turn related to the ef- 
fective co-moving density of Ly-limit absorbers (A x riabs) 
as well as their bias (i?abs) 

CC [nabs(l + ^)'A X Sabs] • (2) 

We can also relate the average co-moving density of sources 
to the collapsed fraction of baryons (Fcoi) in dark-matter 
halos above some minimum mass (that reflects the thresh- 
old for cooling and fragmentation of the gas into stars 
or accreting black holes). We assume that ionizing pho- 
tons were generated by starbursts, and therefore expect 
that the luminosity density of ionizing photons was pro- 
portional to the star-formation rate. We further assume 
that the starbursts formed primarily out of newly collaps- 
ing IGM. As a result the star-formation rate is assumed to 
be proportional to the derivative of the collapsed fraction 
risrc oc (dFcoi/dt), so that 

J.ocl^^^flld + Fc). (3) 

^abs -^abs 

If both sources and absorbers reside in halos of the same 
mass, then Babs = -Ssrc- This may be the case once the gas 



in mini-halos below the cooling and fragmentation thresh- 
old, has been photo-evaporated. However at high redshifts 
just after the overlap of H II regions, the mini-halos might 
provide a significant contribution to the population of Ly- 
limit absorbers (Barkana & Loeb 2002; Iliev, Shapiro, & 
Raga 2005). In this case Bgic > -Babs- These bias factors 
are expected to evolve with time; Bgrc may have a slow 
redshift dependence through a changing minimum mass 
for galaxy formation while i?abs may have a value that 
changes with redshift if the typical absorber mass c;lianges 
with time (for example if mini-halos are being evaporated) 
and so Babs could be a function of the evolving J^. 

Suppose that riabs is some general function of -Fcoi and 
Ji,. Ly-limit systems are expected to shrink in physical 
cross-section as their lower-density envelopes are ionized 
by a growing ionizing background. The effective density 
of absorbers should therefore be a function of J^, in addi- 
tion to Fcoi. Moreover since i^coi and {dFco\/dt) are both 
monotonic there is a unique one-to-one correspondence be- 
tween the values of Fcoi and {dFco\/dt) as a function of 
time. Equation (3) is therefore an implicit equation for 
as a function of Fcoi. The above considerations lead to two 
interesting conclusions. First, we find that J^, should be in- 
dependent of density contrast at a fixed value of collapsed 
fraction; and second, the average ionizing background Ji, 
within a region of radius equal to Rmip should be a function 
only of the value for i^coi- With these points in mind we can 
estimate the scatter in the redshift at which a threshold 
value of Jj/ is reached, and from it find the corresponding 
scatter in the value of Jj, at fixed redshift. 

2.1. Fluctuations due to cosmic variance 

Large-scale inhomogeneity in the cosmic density field 
leads to structure-formation that is enhanced in over-dense 
regions and delayed in under- dense regions. The resulting 
cosmic variance in the redshift at which a critical value of 
i^coi is reached within regions of size R may be calculated 
as 



(4) 



where is the r.m.s. amplitude of the linear density 
field smoothed over spheres of radius R, and Sc is a crit- 
ical overdensity for collapse (oc {1 + z) at high redshift; 
Barkana & Loeb 2004; Bond et al. 1991). Since the ion- 
izing background within a given region depends on the 
value of the collapsed fraction, we can relate {5z1^) 2 to 
the cosmic variance in redshift where a critical value of 
Ju is obtained. The value of Jv is independent of A. We 
therefore use the mcan-free-path averaged over all over- 
densities as the appropriate length scale over which the 
cosmic variance is calculated, i.e. R = i?mfp/2. 

The left panel of Figure 1 shows the scatter in the red- 
shift at which a critical value of J^, is obtained due to cos- 
mic variance, {Sz"^^) (solid lines). We would like to convert 
this scatter in redshift to a scatter in the optical-depth to 
Lya absorption at a given redshift. At a fixed redshift z, 
the local value of Lya optical depth (n) scales as 



(5) 



where pui is the density of neutral hydrogen (Barkana & 
Loeb 2004). The effective optical depth (r) is measured 
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Fig. 1. — Left panel: Contribution to scatter in rcdshift from cosmic variance (5x^^)2 (solid line) and finite light travel time (<52caus)' 
(dotted line) as a function of the mcaii-frcc-path for ionizing photons in co-moving Mpc at a redshift of 5.5. Right panel: The value of 
dlnrigi^s/dz — din {dF^i / dt) / dz as a function of redshift. The dashed, solid and dotted lines correspond to 7 = 2, 3, and 4 respectively. In 
calculating F^^l we have assumed a minimum virial temperature for galaxies of 2 x lO^K. 



over a path-length through the IGM, and may be com- 
puted from T = — logF where F = (e""^') is the trans- 
mission averaged over the path-length. While the ionizing 
background will be sensitive to cosmic variance in density 
contrast calculated in 3-dimensional volumes (A) on the 
scale -Rmfp, the absorption spectra of quasars probe only 
lines-of-sight through these volumes. The effective optical 
depth averaged over a length scale i?mfp therefore depends 
on the square of the density averaged along the line-of- 
sight (Lidz, Oh & Furlanetto 2005). Here we neglect den- 
sity fluctuations on small scales and approximate the av- 
erage of the density contrast squared using the smoothed 
density contrast (A^) calculated in a thin cylinder at red- 
shift z. This procedure underestimates the fluctuations 
in transmission (Lidz et al. 2005). Note that although t 
varies with z due to the growth factor of density perturba- 
tions, we focus here on the scatter in t at a fixed z. The 
fractional change in the value of the optical depth relative 
to the average, given a density contrast (A^), and a delay 
{5z) in the redshift where the critical collapse fraction is 
reached is 

J V (^av) 



5,= 



= A 



- 1 



dz 



dz 



- 1 



(6) 



where Tav is the average value of optical depth at Zav. 
But J,j oc (dFcoi/cJt)n^j;^g(Bsrc/-Babs)(l + ^Q) [equation 3]. 
Within this expression we can only compute the evo- 
lution of collapsed fraction, and therefore assume an 
evolution specified by a parameter 7 and of the form 
^^abs(-Babs/-Bsrc)(l + Fq) cx (1 + z)^ for the remaining 
quantities. Since we are considering a fractional change 
in redshift that is much smaller than unity a power-law 
approximation is an appropriate choice. With this addi- 
tion we get 



1 



dln(rfFcoi/rft) 



7 



1. (7) 



dz (1 + 2:). 

We note that equation (7) is a first order expression 
for This choice may lae justified by the fact that 



dln(dF.3oi/d*) 



has an evolution which is nearly lin- 



ear with redshift (see figure 1). 

At this point we specify the value of 7. One can ob- 
serve the average number of Ly-limit systems per unit 
CO- moving length along the line-of-sight (i?LL)- Given 
the observed distribution of column densities for Ly-limit 
systems, this distance can be related to the mean free 
path for Ly-limit photons -Rmfp = -Rll/a/tt (Miralda- 
Escude 2003). Storrie-Lombardi et al. (1994) measured 
the abundance of Ly-limit systems per unit redshift at 
= 4, dNi,i,/dz ~ 3.3, leading to a co- moving spacing 
of i?LL = [(1 + z)c{d,t/ dz) /Z-i] for Ly-limit systems along 
the line-of-sight. The resulting mean- free-path for Ly-limit 
photons is i?nifp ^ 120 co-moving Mpc at z = 4. This 
value provides an upper limit on the mean-free-path for 
Ly-limit photons at higher redshifts when the cosmic gas 
is denser and less ionized. The evolution of the number 
of Ly-hmit systems per redshift is dN-n^/dz oc (1 -I- z)^ 
with [i ~ 1.5 (Storrie-Lombardi et al. 1994). The evolu- 
tion in the number of lines per unit co-moving length x is 
therefore 



LL 



OC {l + zY 



,dN, 



LL 



cd,t 
dz 



oc(l + z)'^+i, (8) 



dx dz 
where c is the speed of light and we have neglected the in- 
flucmcc of the cosmological constant at high redshift. The 
density of absorbers is proportional to dNi,i^/dx, hence we 
flnd riabs oc {l + z)^~^'2 , and thus 7 = 1. The value of (3 
is meastircd at z = 4. However mini-halos may be impor- 
tant at high redshift so that an extrapolation of absorption 
properties from low to high redshift is not justifled. We 
therefore consider a range of values for /3 which allow for 
a range of evolutionary properties of the absorbers as well 
as their bias, and of the contribution of quasars to the 
ionizing background. In the right panel of Figure 1 we 
plot the values of the difference between dlnriahs/dz and 
d\n{dFcoi/dt)/dz as a function of redshift for (3 = 0.5, 1.5 
and 2.5 (7 = 2, 3 and 4). 

In a region of density contrast A, we can compute the 
typical delay due to cosmic variance {Sz^^)^ (see equa- 
tion 4). However to estimate the typical fluctuations in the 
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effective optical depth we must also compute the variance 
in overdensity among lines of sight through the density 
field of length -Rmfp- We calculate the power-spectrum of 
fluctuations in cylinders of length L and radius R to sec- 
ond order in the wave number k (Kaiser & Peacock 1991; 
Baugh 1996) 

Pm{k) = J^t[ dyyP^(y)e-^'^-^'^^'l\ (9) 

where Pb (fc) is the linear baryonic power-spectrum, which 
may be approximated as Ph{k) — P{k){l + k^Rf}^^ in 
terms of the cold dark matter power-spectrum P{k) and 
the filtering scale Rf for the associated reionization history 
(Gnedin & Hui 1998). The variance in on a scale R 
with wavenumber k^. follows from 

ffcR 

(10) 



1 r"- 

= 27r y dkPid{k). 



The radius of the cylinder R is set by the size of the quasar 
emission region, which is much smaller than the filtering 
scale. Our results are therefore insensitive to the precise 
choice of R. To estimate typical fluctuations in the opti- 
cal depth from cosmic variance we assume and A to 
be approximately uncorrelated and combine the 1-sigma 
fluctuation in delay, {{Sz'^^)^^^ with the 1-sigma fluctua- 
tion in density contrast, (A^ = 1 ± a^^zD, where D is the 
growth factor). The typical positive (-I-) and negative (-) 
fluctuations in optical depth due to cosmic variance may 
therefore be estimated as 

'dln{dFcoi/dt) 



1± 
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1. (11) 



dz (1 + -2) . 

Figure 2 shows the typical positive (+) fluctuation 
((5^ ,,^)2 , which is monotonic in i?mfp. 

2.2. Fluctuations due to the finite light travel time 

Following the appearance of sources with luminosity 
density A x Usrc, the ionizing intensity will grow to its full 
value Jj/ over a timescale i?mfp/[c(l + z)]. The time when 
a critical J,^ is achieved is affected by this light propaga- 
tion delay (with ionizing sources emitting their radiation 
at a retarded time). However, there may also be a scatter 
in the time where the critical is achieved from place 
to place within the region of co- moving size i?mip. This 
scatter will depend on the distribution of sources. For 
example, if the sources were clustered on a scale similar 
to the mean-frec-path, then within a region there will be 
a scatter of -Rmfp [(1 + z)cdt/dz]~''' in redshift where the 
critical is achieved. However if the clustering is on a 
much smaller scale then the scatter will be smaller. 

We can estimate this scatter by considering the clus- 
tering length of sources, Rq (which is ~ 4 co- moving 
Mpc for dark-matter halos with a virial temperature of 
Tvir ~ 2 X lO^K at z ^ 5.5). The photons must travel a dis- 
tance Rmfpi^ + z)~^ before the critical J^, is reached. The 
scatter in arrival times for photons coming from sources 
clustered on a scale Rq at a distance -Rmfp may be esti- 
mated as the crossing time of the smaller chistcring length 
(Ro/c) divided by the square-root of the number (A^ciuster) 
of clustered regions at a distance i?mfp (where the last fac- 
tor corresponds to Poisson fluctuations in the distances 



of the clustered regions). Hence the variance in redshift 
where the critical is reached due to the light propaga- 
tion delay will therefore be 

-1 



i?0 



V^cluster(l + z) 

where A^ciuster is estimated as 



cdt 
dz 



4Tr{Rrnip/2fRo _3 fRn,tp\ 

iVcluster ~ ^^^3^3 ' 4 [ J^^ J 



(12) 



(13) 



The corresponding redshift variance due to the light prop- 
agation delay {Sz'^^^g)^ is shown on the left panel of Fig- 
ure 1. At a fixed density contrast, the finite light travel 
time will introduce a scatter in the value of Sr of 



T.caus/ 



1 



dln{dFcoi/dt) 
dz 
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('^^caus) ■ 



- 1- (14) 



(1 + ^). 

We find that the fluctuations in optical depth due to this 
effect are much smaller than those due to cosmic variance 

(sec Figure 2). 

2.3. Fluctuations due to the finite number of ionizing 

sources 

The Press-Schechter formalism predicts a probability 

distribution (and hence a scatter) for the number of halos 
within a region of a given matter overdensity. In this paper 
we use the derivative of the collapsed fraction of baryons 
(up to a proportionality constant) as a proxy for the star 
formation rate. This quantity depends on the overdensity 
A. As a result there is a scatter in the redshift where 
a given collapsed fraction is reached among representa- 
tive regions of the universe, or equivalently there is scatter 
among the collapsed fractions within different regions at 
fixed time. In this paper we are concerned with regions of 
overdensity A and finite size -Rmfp. Therefore, in addition 
to the fluctuations in the density and ionizing radiation 
introduced by this cosmic variance (and by the light prop- 
agation delay) , we expect a spatial variation of the ioniz- 
ing intensity due to fluctuations in the number of sources 
within the finite region. We assume that these fiuctuations 
may be represented by Poisson noise in the number of ion- 
izing sources within a region of size .Rmfp. We estimate the 
number of galaxies using the Press-Schcchtcr (1974) mass 
function dnps/dM (with the modification of Shcth & Tor- 
men 1999) and evaluate the number per logarithm of mass 
at the minimum mass for galaxy formation (-Mmin)- The 
number of ionizing sources at a particular epoch is given 

by 



dn 



PS 



tit 

tn dliiM 



R 



mfp 



(15) 



where t\t is the source lifetime and t-^ is the Hubble time 
at redshift z. In estimating the fluctuations due to varia- 
tions in the source numbers we consider two limiting cases. 
First, we adopt txt/t-a ~ 0.01 corresponding to starburst 
lifetimes of ^ 10^ years at z ^ 5.5 (the lifetimes of massive 
stars) and a minimum mass corresponding to the Jeans 
mass in a photo-ionized IGM (halo virial temperature of 
~ 2 X lO^K). This case corresponds to the smallest source 
numbers and hence the largest level of Poisson fluctuations 
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R^^p (Mpc CO — moving) 

Fig. 2. Contributions to fluctuations in Sr- The thiclj lines show the fluctuations in optical depth induced through cosmic variance, 
{S'^^^)^^^. The dashed, solid and dotted lines correspond to 7 = 2, 3, and 4 respectively, and for calculation of F^^i we have assumed a 
minimum virial temperature for galaxies of 2 X lO^K. The values of (<5^ Poisson)^^^ ^""^ shown by the light lines; the solid and dashed lines 
correspond to tit/tii = 0.01 and tit/tn = 0.1, combined with minimum virial temperatures for galaxies of 2 X lO^K and lO^K respectively. 
We also show the values of (^rcaus)^^^ (dark lines with dashed, solid and dotted curves corresponding to 7 = 2, 3, and 4 respectively). 
(■^T.caus)^^^ is much smaller than (5^ Poisson)^''^ {^t,cv}^^^- each case we have evaluated the fluctuations at z = 5.5. 



in J^. Second we consider t\t/tu ~ 0.1 corresponding to 

starburst lifetimes of ^ 10^ years (of order the orbital 
time for a circular orbit in a star-forming disk that oc- 
cupies ~ 10% of the virial radius of its host halo). This 
lifetime is combined with a minimum halo mass having 
a virial temperature of ~ lO^K, corresponding to cooling 
threshold of atomic hydrogen (so that the filtering scale is 
assumed not to have had time to respond to the re-heated 
IGM) . This case corresponds to the largest source numbers 
and hence to the smallest fluctuations in J^. The Poisson 
noise in the number of galaxies is AA^s^c = V-^src, so that 
the fluctuations in the optical depth induced by Poisson 
fluctuations in the number of sources is 

where we have used the relation Ji, cx Nsrc, and evaluated 
the fluctuation range of r corresponding to the 1-sigma 
level of fluctuations in the number counts of sources. The 
resulting curves are plotted as a function of i?mfp in Fig- 
ure 2. 

2.4. Fluctuations in the optical depth 

The contributions to the redshift variance introduced by 

causality, by cosmic variance and by source count fluctu- 
ations are independent of each other. The positive (-I-) 
and negative (-) ranges of allowed Sr may therefore be 
estimated by adding these contributions in quadrature, 

iS'S' = ±v/(^?,cv) + (^?,caus) + ('^?,Poisson)- (17) 



3. RESULTS AND DISCUSSION 

The upper and lower limits that bracket the typical 
range of Sr at z = 5.5 are plotted in Figure 3. The dashed, 
solid and dotted lines correspond to 7 = 2, 3 and 4 re- 
spectively. The left and right panels show the cases of 
tit /in = 0.01 and 0.1, combined with minimum galaxy 
masses corresponding to virial temperatures of 2 x lO^K 
and lO'^K respectively. Figure 3 shows that cosmic vari- 
ance introduces fluctuations in r that are of order unity 
over a range of length scales. As the mean free path in- 
creases the fluctuations become smaller in line with the 
smaller cosmic variance on large scales. The fluctuations 
in introduced by fluctuations in the number of sources 
provides a small contribution. This is in contrast to the 
situation for double reionization of helium at lower red- 
shifts where the ionizing sources (quasars) reside in rarer 
more massive systems (Bolton et al. 2005). 

The finding that cosmic variance should result in signif- 
icant fluctuations in r is qualitatively insensitive to varia- 
tions in the value of 7 (which parameterizes unknowns re- 
garding the absorber properties) over a wide range. This 
insensitivity arises partly because the collapsed fraction 
evolves much more rapidly than a power-law with redshift 
(because the minimum galaxy mass lies on the exponen- 
tial tail of the Gaussian density distribution), and partly 
because the fluctuations are dominated by the term in 
the evaluation of r. 

We may take the upper limit for Rmfp after reionization 
to be 120 co-moving Mpc, which is the value measured at 
2; ~ 4 (Miralda-Escude et al. 2003), while the lower limit is 
set by the size of bubbles at the end of the overlap epoch 
(~ 70 co-moving Mpc at z ^ 6 and larger if reionization 
completed earlier; Wyithe & Loeb 2004b). The shaded 
region in Figure 3 delineates the range for different values 
of Rmfp assuming 7 = 3. We find our main conclusion. 
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Fig. 3. — The upper and lower values of the 1-sigma range for fluctuations in St- The dashed, solid and dotted lines correspond to 7 = 2, 3 
and 4 respectively. The grey portion corresponds to the range spanned by different values of iimfp given 7 = 3. The grey region is bounded 
from the left by the boundary at 70 co-moving Mpc corresponding to the bubble size at the end of then reionization, and from the right by 
the ionizing photon mean-free-path at z ~ 4. We show cases with z = 5.5, and take populations of ionizing sources having tit/tn = 0-01 (left) 
and 0.1 (right), combined with minimum galaxy masses corresponding to virial temperatures of ~ 2 x 10* K and ~ 10* K respectively. The 
values of 5r at z ~ 5.6 measured along the lines-of-sight to five high-redshift quasars (White et al. 2003) are over-plotted. The horizontal 
error-bars on these observations indicate the possible range of iimfp . 



that cosmic variance introduces large fluctuations in the 
transmission of the IGM at z ~ 5.5, holds over this full 
allowed range of Ly-limit mean- free-paths. The reioniza- 
tion rcdshift itself is also inhomogcncous, as evidenced by 
the transmission of flux in the Lya troughs of some but 
not all of the z > 6 quasars (White et al. 2003; Oh & 
Furlanetto 2004). Our model becomes applicable at the 
redshift where most lines of sight probe the post-overlap 
IGM. Finally, we note that our results do not depend on 
the reionization redshift having occurred near z 6 but 
apply to earlier reionization epochs as well. 

3.1. Scatter in the transmitted intensity 

The fluctuation in the relative intensity (/ a e~'^) trans- 
mitted through the IGM (Sj = (I — lav) /lav) is related to 
the fluctuation in r through 



lav e 
= e-''--^- - 1. (18) 

Note that in regions where the optical depth is reduced 
{Sr < 0) the exponential is positive and the value of Si 
can become very largo. For values of optical depth larger 
than unity the exponential acts as a non-linear amplifier 
of any inhomogeneities in the IGM which skews the distri- 
bution of transmitted intensity. Cosmic variance therefore 
leads to a highly asymmetric distribution of transmitted 
flux, even if the distribution of Sr is a symmetric Gaussian. 
For large values of r ~ 3-5 immediately following reion- 
ization, the variations in the transmission span an order of 
magnitude. At later times when the optical depth becomes 
small, equation (18) becomes 

Si ~ -TaySr. (19) 

Thus at small values of optical depth the distribution of 
transmitted intensity becomes symmetric. 

3.2. Comparison with observations 



The results presented in Figure 3 may be compared with 
existing observations of high redshift quasars. Becker et 
al. (2001) and White et al. (2003) calculated the effective 
Lya optical depth in bins of width Az = 0.2 along the 
lines of sight to five high redshift quasars. The resulting 
values of Sr for these five lines of sight in the redshift bin 
5.5 < z < 5.7 are plotted over the grey region in Fig- 
ure 3 for comparison. At 2: ~ 5.6 White et al. (2003) find 
Tav ~ 2.5 and fractional fluctuations in optical depth with 
values up to (r — Tav)/Tav ~ 0.25. Since there are only 5 
lines of sight we would expect this observed range to be 
roughly equivalent to the 1-sigma fluctuations plotted in 
figure 3, for which we find ((5^) = ~ 0.1 — 0.4 depending on 
the value of the mean- free-path and the nature of the ioniz- 
ing sources. Inspection of Figure 3 therefore suggests that 
the model is in good agreement with with current data. 
Djorgovski (2004) has also presented the fluctuations in 
optical depth along five lines of sight within Gaussian red- 
shift windows of variance = 0.05. This smaller redshift 
window yields fractional fluctuations in optical depth with 
values up to (r — Tav) /Tav ~0.5. Atz~5.5, the redshift 
bin widths of Az = 0.1 — 0.2 correspond to co-moving line- 
of-sight distances Raz = (1 + z)cdt/dzAz ~ 50 — lOOMpc. 
Thus the observations sample the IGM over a length-scale 
comparable to i?mfp, and so are fortuitously matched to 
our calculation. Of course at redshifts where i?mfp < 
50 — lOOMpc, our estimate of (^^) 2 may be compared with 
observation by averaging the predicted fluctuations over a 
number ~ [(50 — 100Mpc)/_Rnifp] of regions. 

3.3. Relative absorption of Lya and Ly(3 photons 

The relative level of absorption for Lya and Ly/3 (or 
Ly7 etc) photons is dictated by the ratio of the product 
of oscillator strengths and resonance wavelengths for the 
two transitions. Assuming a homogeneous IGM this re- 
sults in a Lya optical depth that is 6.25 times the Ly/3 
optical depth (Songaila 2004). However the ratio is de- 
creased in an inhomogeneous IGM because on small scales 
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transmission is dominated by low-density regions, while 
the transition with the lower cross-section for absorption 
is able to sample a larger range of densities. Oh & Furlan- 
etto (2004) find that values of 2-2.5 are more realistic for 
the ratio of Lya to Ly/3 optical depth. 

The constant of proportionality relating r and Ji, will 
therefore be dependent on the oscillator strength and res- 
onance frequencies. This constant cancels in calculation of 
(5t, and so the range of fractional variation in optical depth 
will bo the same for absorption in Lya and Ly/3. On the 
other hand the range of optical depth values is given by 
Tav(<5^)^. As a result, the absolute range of optical depths 
observed for Lya absorption will be larger than the range 
observed for Ly/? absorption by a factor of ^ 2 — 6. Sim- 
ilarly, while the relative fluctuations in the optical depth 
are the same for all transitions, we flnd that the relative 
fluctuations in observed intensity 5i will be much larger 
for Lya than Ly/3 photons due to the non-linearity of the 
exponential relation between optical depth and transmis- 
sion. In particular, if the Lya optical depth is smaller than 
unity, then Sj ^ —Ti^vSt and the relative fluctuations in in- 
tensity may be calculated directly from the ratios of the 
optical depths among the different transitions. In cases 
where the Lya optical depth is larger than unity, the vari- 
ation in 6i among the different transitions is larger due to 
the presence of the exponential in equation (18). 

4. SUMMARY 

In this paper we have calculated the magnitude of the 
large scale fluctuations in the optical depth to absorption 
of Ly-series photons that are introduced through variance 
in the cosmic density field. We find that cosmic variance 
introduces fractional fluctuations in the optical depth that 
are of order unity for mean-free paths of ionizing pho- 



tons i?nifp ^ 150 co-moving Mpc. This finding explains 
the variation seen in the Lya optical depth along different 
lines of sight through the IGM at ^; < 6 (Becker et al. 2001; 
White et al. 2003). 

Our model was not constructed to predict the value of 
optical depth or to associate a value of the optical depth 
with the mean- free-path. Rather, it predicts the trend 
that the relative fluctuations in optical depth will drop as 
the ionizing photon mean-frec-path grows. The mean-free- 
path increases with cosmic time, and so the range of opti- 
cal depth will therefore decrease towards low redshift. In 
addition, while our model predicts a distribution of optical 
depths that is roughly symmetric about the average, the 
exponential dependence of transmission on optical depth 
implies a distribution of transmission that will be highly 
asymmetric. This asymmetry will become smaller as the 
average value of optical depth declines at late times. Fi- 
nally, our model predicts that the absolute fluctuations in 
the Lya optical depth should be greater than in the Ly/3 
optical depth. This translates to relative fluctuations in 
the transmission through the IGM that are much larger 
for Lya than for Ly/3. 

As larger numbers of high redshift quasars arc discov- 
ered through surveys like the Sloan Digital Sky Survey, the 
statistics of absorption averaged over large scales will be 
measured along many lines of sight. The results of this 
paper suggest that the effects of cosmic variance on the 
post-reionization epoch will be directly observed in these 
studies. 

The authors would like to thank Xiaohui Fan for point- 
ing out an error in an early version of this manuscript. 
JSBW acknowledges the support of the Australian re- 
search council. This work was supported in part by NASA 
grants NAG 5-13292 and NNG05GH54G for AL. 



REFERENCES 



Barkana, R., Loeb, A., 2002, ApJ, 578, 1 

. 2004. ApJ, 601, 64 

Baugh, C. M. 1996, MNRAS, 282, 1413 
Becker, R. H., et al. 2001, AJ, 122, 2850 

Bolton, J. S., Haehnelt, M. G., Viel, M., & Carswell, R. F. 2005, 

MNRAS, submitted; astro-ph/0508201 
Bond, J.R., Cole, S., Efstathiou, G. & Kaiser, N., 1991, ApJ, 379, 

440 

Ciien, H.-W., Prochaska, J. X., Bloom, J. S., & Thompson, I. B. 

2005, ApJL, submitted; astro-ph/0508270 
Djorgovski, S.G., 2004, in Proc. Tenth Marcel Grossmann Meeting 

on General Relativity, cds. M. Novello, S. Ferez-Bergliaffa and R. 

Ruffini, astro-ph/0469378 
Fan, X., Narayanan, V. K., Strauss, M. A., White, R. L., Becker, 

R. H., Pentericci, L., & Rix, H.-W. 2002, A J, 123, 1247 
Fan, X., et al. 2004, AJ, 128, 515 

Furlanctto, S. R., Hernquist, L., & Zaldarriaga, M. 2004, MNRAS, 

354, 695 

Gnedin, N. Y., & Hui, L. 1998, MNRAS, 296, 44 

Iliev, I. T., Shapiro, P. R., & Raga, A. C. 2005, MNRAS, 361, 405 



Iliev, I.T., McUcma, G., Pen, U.-L., Merz, H., Shapiro, P.R., Alvarez, 

M., 2005, astro-ph/0512187 
Kaiser, N.. & Peacock, J. A. 1991, ApJ, 379, 482 
Kohler, K., Gnedin, N.Y., Hamilton, A., 2005, astro-ph/0511627 
Kohler, K., Gnedin, N. Y., Miralda-Escude, J., & Shaver, P. A. 2005, 

ApJ, 633, 552 

Lidz, A., Oh, S.P., Furlanetto, S., 2005, astro-ph/0512427 

Mcsingcr, A., & Haiman, Z., 2004, ApJL, in press, astro-ph/0406188 

Miralda-Escude, J. 2003, ApJ, 597, 66 

Oh, S. P., & Furlanctto, S. R. 2005, ApJ, 620, L9 

Press, W., Schcchtcr, P., 1974, ApJ., 187, 425 

Sheth, R. K. & Tormcn, G. 1999, MNRAS, 308, 119 

Songaila, A. 2004, AJ, 127, 2598 

Spergel, D. N, et al., 2003, Astron. J. Supp., 148, 175 
Storrie-Lombardi, L. J., McMahon, R. G., Irwin, M. J., & Hazard, 

C. 1994, ApJ, 427, L13 
White, R., Becker, R., Fan, X., Strauss, M., 2003, Astron J., 126, 1 
Wyithe, J.S.B., Loeb, A., 2004a, Nature, 427, 815 
Wyithe, J. S. B., Loeb, A., 2004b, Nature, 432, 194 



